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PREDICTION OF SOLAR PROTON DOSES IN THE APOLLO 


COMMAND MODULE AND LEM 
Ely Terence M. Vinson and Manuel D. Lopez 


ABSTRACT 


A method of predicting the solar proton dose inside the Apollo 
Command Module and LEM for various stay times in each vehicle is dis- 
cussed. Using analytical expressions to obtain the time dependent 
flux for 5^ recorded proton events, a computer program generates pro- 
bability distributions of dose for various missions. From these dis- 
tributions, curves of dose as a function of stay time for each vehicle 
are constructed. 


INTRODUCTION 


In a previous NASA TN by Modisette, Vinson, and Hardy (ref. l), 
a method of determining the solar proton environment for missions of 
from one week to two years was discussed. In this method, the time- 
integrated fluxes of the 5^ largest recorded events in the upper half 
of the last solar cycle were used as input data. A computer program 
was devised which calculated the total proton flux encountered during 
all possible missions of the desired length. The proton events were 
treated as ’’spikes" of flux; that is, the time-integrated flux for an 
event was assumed to be encountered in a single day., The total en- 
countered fluxes were assumed to fit the normal probability distribu- 
tion. Probability distributions of encountered flux were made up for 
each mission length, and from these, curves of flux as a function of 
mission length at certain probability levels were obtained. 

This approach is valid for missions of the order of a week or 
longer, but for missions which are shorter than the actual duration 
of an event, a more sophisticated technique is necessary. 


SOLAR PROTON DOSE ANALYSIS 
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General Method 

The method described in reference 1 can be modified so that it can 
be used to predict doses for ^missions”, or stay times , which are of the 
order of a few hours long. Instead of considering every day in the six- 
year period covered by the data as the beginning of a mission, every 
hour of each day is so considered. A computer program then calculates 
the total flux encountered during each mission and the flux is converted 
to dose. The program then generates a probability distribution of doses 
encountered for that particular stay time. 


Input Data 


It was necessary in this analysis to describe as accurately as pos- 
sible the time history of each proton event considered. To facilitate 
calculating hourly doses from the proton events, analytical expressions 
for proton flux were sought which reasonably filled available data. The 
rise time, event duration, peak and time- integrated intensities, as de- 
termined analytically, were required to match those of actual events at 
a rigidity of .2391^V (30 MeV). 

The expression used to determine particle flux as a function of 
time t and rigidity P is: 


n 


-1 


[• 


F(P,t) - At** exp I - BPt n - Ct 


-n] 


J 


a) 


, are determined from the analytical fit 

n 


The constants. A, B, C, n 
for each proton event. The time -integrated flux, J F(P,t)dt, is pro- 
portional to e " < V^ r where d is a constant, and in most cases a value 
for d can be found such that there is a favorable agreement with the 
time -integrated proton fluxes as given by Webber (ref. 2). 


A computer program was written which evaluates the constants A, B, 
C, and n from proton event data. F(P,t) was then evaluated on an 
hourly basis giving the characteristic spectrum slope P_ and the in- 


,, ° st, 

tensities for rigidities greater than 0, .1 366, .2391? and .4448"5£8f* * 

The hourly P and the flux for P greater than . 2391 BeV' were then 

o py 

used to calculate the corresponding dose rates. The dose rate calcula- 
tions used normalized proton dose versus P relationships for the 
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command module and LEW as calculated by Hardy (ref, 3). For practical 
purposes, the dose rate calculated for a particular hour was considered 
to be the total dose for that hour. 


Analysis of Results 

The object of this study was to obtain probability distributions of 
the doses encountered during various stay times in the command module and 
the LEM, From these distributions, curves of dose versus stay time at 
various probability levels can be constructed. The distribution of dose 
received in the LEM is of greater interest because the LEM is a lightly 
shielded vehicle and it is to be expected that most of the dose will be 
received either in the LEM or on the lunar surface in the space suit. 

The probability that a solar proton dose will be greater than D 
rads is given by the following equation which assumes that the doses are 
normally distributed: 


where 



( 2 ) 


« mean of dose distribution for mission length L 

* standard deviation of dose distribution for mission length L 

For a given mission length, the probability P g of encountering a dose 
of any size is given by: 

♦ 

P a missions which encounter a dose /,% 

e * total number missions 

The probability of encountering a dose greater than D is therefore 
given by: 

P D « P e . P(z) (4) 

P is a constant for a given mission length. An arbitrary value such 
e 

as ,01, .005, or .001 can be assigned to P^, so that equation (4) can 
be solved for P(z) and hence D can be obtained. 
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Figure I and 2 show the probability distributions of skin dose 
for stay times in the LEM of 24 and 48 hours. Figures 5, 4, and 5 cho" 
the distributions for stay times in the command module from 48 hours 
to two weeks. The straight lines in the figures were dravn by means of 
a normal fit using the mean and standard deviation of each set of data 
points. The data appear to fit a normal distribution reasonably well, 
except for the longer missions in the command module. This may be due 
in part to the manner in which the data were grouped. As in refer- 
ence (l), the doses obtained were grouped according to the length of 
the mission in order to avoid redundancy in the large number of data 
points obtained. Actually, the averages of the groups of doses rather 
than the doses themselves are plotted. This, of course, tends to make 
the distribution more nearly normal. 

Figures 6 and 7 show skin dose as function of time in the command 
module and LEM for three probability levels. From these curves one can 
predict the dose for a combination of time spent in both vehicles. For 
example, suppose an astronaut spent 72 hours in the command module, then 
24 hours in the LEM, then 72 hours in the command module. For each 
stay time, what would the maximum dose be with probability ,01? From 
figures 6 and 7 the following numbers are obtained: 

72 hours in the command module . l4 rads 
24 hours in the LEM - 163 rads 

72 hours in the command module l4 rads 

193 rads 

For each of the above stay times, the probability Is 1,0 - .01 * .99 
that the given dose will not be exceeded. The probability of not ex- 
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ceeding the dose during all three times taken in sequence is , 99 * .97* 

Therefore, the probability of exceeding 193 rads for this sequence of 
stay times is 1.00 - .97 * , 03 , 

This analysis is subject to possible error from three important 
sources. These are (l) the dose calculations for the two spacecraft, 

(2) the expressions used to generate the time -dependent fluxes for 
each event, and ( 3 ) the assumption that the doses are normally dis- 
tributed. The second source of error is probably the most serious. 
Although no claim is made for the accuracy of the expressions used, 
they are probably the best that could be obtained, given the limited 
amount of data available on the fluxes and spectra of the actual events. 


SUMMARY AND CONCLUSIONS 

The following table shows a comparison of some of the results ob- 
tained in this study to some results obtained by an alternate method. 



Command module 
stay time, hours 


48 

100 

l40 


LEM stay 
time, hours 

6 


25 


50 

60 


TO 


? TWO METHODS OF DOSE PREDICTION 


.01 probable dose, rads 

Factor 

Method A 

Method B 

difference 

51.5 

7-8 

4.03 

42 

21 

2.00 

44 

55 

1.55 

.01 probable dose, rads 

Factor 

Method A 

Method B 

difference 

9 

47 

5.22 

190 

170 

1.12 

450 

518 

1.55 

500 

574 

1.54 

550 

450 

1.28 
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In method A, the Integrated flux of a .01 probable event as determined 
in reference 1 was filled to a model time-variant spectrum suggested by 
Bailey (ref. 4), The dose calculations were performed by Hardy (ref, $). 
Under method B are listed ,01 probable doses as calculated by the method 
of this paper. 

The doses obtained by method B are consistently lower than those 
obtained by method A with the exception of the six-hour dose in the LEM. 
The reason method A gives a lower dose in this case is probably because 
this method calculates the accumulated dose as a function of time after 
onset of the event, and in six hours the event has not reached its peak 
intensity. Method B, on the other hand, gives the .01 probable dose as 
obtained from the distribution of all six-hour doses, regardless of what 
point in the duration of an event they may occur. 

In general, method A would be expected to give higher doses except 
for very short times, and the data in the table support this. The rea- 
son for this is that method A is based on a model event which has an 
integrated flux greater than any actual event ever recorded. This event 
had an arbitrary duration of 150 hours, which would cause the flux at 
any particular time to be unusually high. Method B uses models of actual 
events, and in most of these events the duration is 225 hours. In most 
cases, therefore, the Intensities will be lower than those used in 
method A. Considering these factors, the agreement between the two 
methods is not bad. One may surmise that if the duration of the model 
event of method A had been longer, the agreement with method B would be 
better. 

At present, this analysis is limited to dosage from solar protons. 

It does not include the effects of alpha particles, which are believed 
to be present in large events in numbers approximately equal to those 
of the protons. lhe input data can be changed to include alpha doses 
without too much difficulty. The dose will then be calculated in REM 
to take into account the different relative biological effectiveness of 
protons and alphas. 
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